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ABSTRACT 


This  report  discusses  the  techniques  and  analysis  for  the  precise 
reduction  of  the  topocentric  and  apparent  places  of  stars  and  artificial 
satellites.  Included  is  a full  discussion  of  the  method  of  star  constants, 
independent  day  numbers,  geocentric  parallax,  parallactic  refraction, 
and  errors.  Also  described,  in  detail,  is  the  analysis  necessary  for 
the  real  time  modeling  of  the  telescope-camera  system,  its  theoretical 
basis,  and  differential  reduction  procedures.  In  addition,  an  original 
method  for  computing  the  distance  of  an  artificial  satellite  from  two 
measurements  of  position  and  one  of  angular  velocity  is  developed. 


ill 


TABLE  OF  CONTENTS 


Abstract  iii 

I.  Introduction  1 

II.  Mean  Place  to  Apparent  Place  to  Topocentric  Place  Reductions  2 

III.  Differential  Reductions  15 

A.  Star  Constants  and  Besselian  Day  Numbers  16 

B.  Independent  Day  Numbers  18 

C.  The  Polar  Regions  19 

D.  Apparent  Place  to  Topocentric  Place  20 

IV.  Ideal  Astronomical  Photography  22 

A.  Standard  Coordinates  22 

B.  The  Method  of  Dependences  29 

V.  Real  Astronomical  Photography  35 

VI.  The  Telescope  Model  41 

VII.  Artificial  Satellite  Reductions  47 

A.  Distance  Estimation  47 

B.  Refraction  52 

C.  Planetary  Aberration  56 

D.  Diurnal  Aberration  56 

E.  Geocentric  Parallax  57 

References  58 


V 


I. 


INTRODUCTION 


The  purpose  of  this  report  is  three-fold.  One  aim  was  the  complete 
elucidation,  through  terms  of  the  second-order  in  small  quantities,  for 
reducing  the  mean  place  of  a star  to  its  apparent  place  and  thence  to 
its  topocentric  place.  Another  reason  was  to  develop  the  basis  for  the 
telescope  modeling  used  in  the  GEODSS  Local  Astrometric  Calibration 
Procedure  and  prepare  the  way  for  the  utilization  of  charge  coupled 
devices  or  charge  injected  devices.  These  two  rely  on  classical  photographic 

astrometry.  In  particular,  the  plate  modeling  used  in  both  the  method 

1 2 
of  dependences  (Schlesinger  ) and  the  plate  overlap  technique  (Eichorn  ) 

are  extended  to  the  real-time  problem  posed  by  artificial  satellite 

reductions.  Lastly,  the  analysis  necessary  for  the  recovery  of  an 

artificial  satellite's  geocentric  position  from  its  measured,  topocentric 

position  is  presented. 

For  the  most  part  the  material  contained  herein  is  a review  and 

summary  of  the  relevant  astronomical  literature.  Original  contributions 

are  scattered  throughout  (principally  Eqs.  (12),  § IIIB,  § VI,  § VIIA) . 

The  epoch-to-epoch  reductions  of  a star's  mean  place  and  its  cataloged 
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mean  place  have  already  been  completely  dealt  with  (Taff  ’ ).  Hence, 

the  initial  data  are  position,  proper  motion,  and  annual  parallax  at  the 
beginning  of  the  nearest  Besselian  solar  year. 
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II.  MEAN  PLACE  TO  APPARENT  PLACE  TO  TOPOCENTRIC  PLACE  REDUCTIONS 


The  reduction  procedures  now  used  for  the  mean  place  to  apparent 
place  reductions  were  developed  by  Bessel.  They  represent  approximations 
of  the  rigorous  formulas  which  retain  all  second  order  terms.  Their 
accuracy  is  further  enhanced  by  restricting  their  applicability  to  time 
intervals  of  ±0.5  year  duration. 

A star's  position  in  the  equatorial  coordinate  system  is  specified 
by  its  right  ascension  and  declination  (i.e.,  longitude  and  latitude). 

The  equatorial  coordinate  system  is  specified  by  its  pole  (the  North 
Celestial  Pole  which  is  the  point  on  the  celestial  sphere  penetrated  by 
the  continuation  of  the  Earth's  axis  through  the  Earth's  north  pole), 
its  fundamental  circle  (the  celestial  equator  which  is  the  great  circle 
on  the  celestial  sphere  described  by  the  projection  of  the  Earth's 
equator  as  seen  from  the  center  of  the  Earth),  its  zero  vertical  (half 
the  great  circle  through  the  poles  of  the  celestial  equator  whose 
intersection  with  the  celestial  equator  marks  the  zero  of  the  longitudinal 
coordinate,  i.e.,  the  vernal  equinox  or  the  First  Point  of  Aries),  its 
heliticity  (right-handed),  and- its  epoch.  The  epoch  is  denoted  by  a 
phrase  such  as  "for  equator  and  equinox  of  1977.0,"  etc.  One  must 
prescribe  the  epoch  for  both  the  celestial  equator  and  the  vernal  equinox 
because  the  equatorial  coordinate  system  is  not  fixed  relative  to  an 
inertial  one. 

The  mean  place  of  a star  is  denoted  by  (o^^j  5^)  • (The  words  place 
and  position  are  used  interchangeably  herein) . For  historical  reasons 
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the  values  of  right  ascension  and  declination  found  in  catalogs  (catalog 

Z, 

mean  place)  include  the  small  effects  of  elliptic  aberration  (Taff  ) . 

As  the  mean  place  of  a star  defines  that  direction  in  which  an  observer 
situated  at  the  solar  system  barycenter  would  view  the  star, 

catalog  mean  place  = mean  place  + e-terms. 

The  true  place  of  a star  is  its  solar  system  barycentric  position  referred 
to  the  true  (i.e.,  actual  at  date)  equator  and  equinox.  It  differs  from 
the  mean  place  by  the  effects  of  precession,  proper  motion,  and  nutation. 
The  apparent  place  of  a star  is  its  geocentric  position  referred  to  the 
true  equator  and  equinox.  It  differs  from  the  true  place  by  the  effects 
of  annual  aberration  and  annual  parallax.  The  topocentric  place  of  a 
star  is  its  position  as  determined  by  an  ideal  telescope.  It  differs 
from  the  apparent  place  by  the  effects  of  astronomical  refraction, 
diurnal  aberration,  and  geocentric  parallax. 

The  two  common  methods  for  performing  the  mean  place  to  apparent 
place  reduction  use  either  (i)  star  constants  in  right  ascension  (a,  b, 
c,  d) , star  constants  in  declination  (a',  b' , c' , d'),  and  Besselian  day 
numbers  (A,  B,  C,  D,  E,  J,  and  J'),  or  (ii)  independent  day  numbers  (f, 
g,  G,  h,  H,  and  i) . C and  D are  sometimes  referred  to  as  aberrational 
day  numbers  since  they  are  used  to  correct  for  annual  aberration  while 
A,  B,  and  E are  used  to  correct  for  precession  and  nutation.  In  addition, 
because  they  represent  terms  of  the  second  order,  J and  J'  are  sometimes 
referred  to  as  second-order  day  numbers.  All  of  these  quantities  are 
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tabulated  in  the  American  Ephemeris  and  Nautical  Almanac,  In  the  following 
A-E,  J,  J',  f,  g,  h,  and  i are  in  seconds  of  arc  and  G and  H (as  well  as 
right  ascension  and  declination)  are  in  radians.  The  proper  motion  in 
right  ascension  (p)  and  the  proper  motion  in  declination  (p ' ) are  in 
"/yr. 

When  all  terms  of  the  second  order  are  included  the  reduction 
procedures  are  of  comparable  accuracy.  However,  the  method  employing 
star  constants  and  Besselian  day  numbers  is  to  be  preferred  when  simultaneously 
reducing  the  positions  of  several  stars  to  the  same  instant  of  time 
while  the  method  employing  independent  day  numbers  is  to  be  preferred 
when  reducing  the  position  of  a single  star  for  several  different  times. 
Moreover,  the  latter  method  is  more  expeditious  when  simultaneously 
reducing  the  positions  of  several  stars  close  together  on  the  celestial 
sphere  to  the  same  instant  of  time. 

If  (a,  6)  denotes  the  apparent  place  and  t is  the  fraction  of  the 
year  to  (from)  the  beginning  of  the  nearest  Besselian  solar  year,  then 

2 

a = a + sinl"[Tp  + Aa  + Bb  + Cc  + Dd  + E + Jtan  6 
o o 

+ (tt/k)  (Cdsece  - Dccose)],  (la) 

= a + sinl"[Tp  + f + gsin(G  + a )tan6 
o o o 

+ hsin(H  + a )sec6  + Jtan^6  ],  (2a) 

o o o 
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(lb) 


6=6  + sinl"[xy'  + Aa'  + Bb'  + Cc'  + Dd'  + J'tan6 

o o 

+ (tt/k)  (Cd ' sec£  - Dc'cose)], 

=6  + sinl"[xy'  + gcos(G  + a ) + hcos(H  + a )sin6 

o o o o 

+ icos6  + J'tan6  ],  (2b) 

o o 


where 


= m/n  + sina  tan6  , 
o o 

a' 

= cosa  , 
o 

(3a) 

= cosa  tan6  , 
o o 

b' 

= -sina  , 
o 

(3b) 

= cosa  sec6  , 
o o 

c' 

= tanecos6  - sina  sin6  , 
o o o 

(3c) 

= sina  sec6  , 
o o 

d' 

= cosa  sin6  , 
o o 

(3d) 

and  Ti  (the  annual  parallax  in  seconds  of  arc)  has  been  set  equal  to  zero 
(for  simplicity)  in  Eqs.  (2).  The  auxiliary  quantities  are  < = the 
constant  of  aberration  (=  20V4958) , m = the  centennial  precession  in 
right  ascension,  n = the  centennial  precession  in  declination,  and  e = 
the  true  obliquity  of  the  eliptic.  If  T denotes  the  number  of  tropical 
centuries  since  1900.0  then. 


30?!2337  + 

0!i8630T  + 8®0  X 10"S^, 

(4a) 

2004V685  - 

0V8533T  - 3V7  x 10"'^T^, 

(4b) 

23°27'8V26 

- 46V845T  - 0V0059T^  + 0'.'00181T^ . 

(4c) 
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The  small  difference  between  the  length  of  a tropical  year  and  a Besselian 

g 

solar  year  (the  latter  is  shorter  by  0.148T)  is  negligible. 

The  independent  day  numbers  are  related  to  the  Besselian  day  number 

via 


f = (m/n)A  + E,  i = Ctane, 

(5a) 

gsinG  = B,  gcosG  = A, 

(5b) 

hsinH  = C,  hcosH  = D. 

(5c) 

The  second-order  day  numbers  are  calculated  (using  the  ± sign  when 
6o  ^ 0)  by 

J = sinl"[(A  ± D)sina  + (B  ± C)cosa  ] [ (A  ± D)cosa 

o o o 


- (B  ± Osina  ] , (6a) 

o 

= sinl"[gsin(G  + a ) ± hsin(H  + a )][gcos(G  + a ) 

o o o 

± hcos(H  4-  a ) ] , ‘ (7a) 

o 

J'  = -sinl"[(A  ± D)sina  + (B  ± C)cosa  ]^/2,  (6b) 

o o 

= -sinl" [gsin(G  + a ) ± hsin(H  + a )]^/2.  (7b) 

o o 


Finally  the  time  of  the  start  of  the  Besselian  solar  year  (e.g.,  the 
instant  when  the  right  ascension  of  the  fictitious  mean  sun  is  18^40™) 
can  be  computed  from  the  expression 
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Jan  0^813516  + 1(24^219878  - ofo00308T)  - [25T]  + LPYR,  (4d) 

where  LPYR  = 1 if  25T  - [25T]  = 0,  0 otherwise.  Thus,  1977.0  = Jan 
0^4626.  This  formula  is  valid  if  T e [0.25,  0.99]. 

As  mentioned  above  these  quantities  are  tabulated  in  the  American 
Ephemeris  and  Nautical  Almanac.  The  dependent  variable  is  (necessarily) 
ephemeris  time  but  tables  in  which  the  mean  sidereal  time  at  Greenwich 
is  the  argument  are  also  included.  Since  the  future  relationship  between 
universal  time  and  ephemeris  time  can  only  be  deduced  from  observations 
not  yet  made,  the  latter  tables  are  not  exact.  The  error  induced  can't 
exceed  iO'.'Ol  unless  |6  I > 84°. 

In  order  to  calculate  the  day  numbers  one  needs  to  know  the  time. 

The  development  of  accurate  clocks,  telescopes,  and  sophisticated  reduction 
techniques  has  led  to  the  discovery  of  irregularities  in  the  Earth's 
rotation.  Hence,  there  are  three  different  universal  times.  UTO  is 
universal  time  as  deduced  directly  from  observations  of  the  stars  and 
the  fixed,  numerical  relationship  between  an  interval  of  universal  time 
and  the  corresponding  interval  of  sidereal  time  (see  below) . UTl  is  UTO 
corrected  for  the  motion  of  the  poles  (also  called  polar  wandering) . 

UTl  represents  the  true  angular  rotation  of  the  earth  and  is  independent 
of  the  observer's  location.  UT2  is  UTl  corrected  for  the  average  seasonal 
variations  in  the  Earth's  rotation  rate  (due  to  polar  cap  melting, 
etc.).  However,  UT2  has  not  been  freed  of  secular  (i.e.,  tidal  friction) 
or  other  Irregular  terms.  Both  the  U.  S.  Naval  Observatory  and  the 
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National  Bureau  of  Standards  have  atomic  clocks  designed  to  reproduce 
UT2.  The  time  they  do  produce  is  called  coordinated  universal  time 
(UTC)  and  is  distributed  by  radio  station  WWV.  Announced  frequency  and 
phase  offsets  keep  UTC  within  ±0?1  of  UT2. 

The  measures  of  mean  universal  and  mean  sidereal  times  are  related 

via 


1 mean  sidereal  day  = 23^56"*4f09054  of  mean  solar  time, 

1 mean  solar  day  = 24^3™56?55536  of  mean  sidereal  time. 

Hence,  when  both  are  in  the  same  units 

1 mean  sidereal  day/1  mean  solar  day  = 0.9972695664 

(=  1/1.0027379093).  When  the  zero  point  is  fixed  (usually  at  O^UT  of 

s 

Jan  0)  these  scale  factors  allow  conversion  with  an  accuracy  of  ±0.01. 

If  more  precision  is  required  or  an  apparent  time  is  needed  the  equation 
of  the  equinoxes  (=  the  difference  between  apparent  and  mean  sidereal 

s 

time;  it's  always  within  ±1  of  0)  must  be  considered. 

Let  (Nj^,  T^)  represent  the  time  of  the  observation.  Here  ^ 0 

is  the  number  of  whole  days  elapsed  since  Jan  1.0  = O^UT  Jan  1.  T is 

u 

UTC,  T^e[0,  24^).  Then  if  A is  the  east  longitude  of  the  observer  (in 

hrs)  and  T (1.0)  is  the  mean  sidereal  time  at  O^UT  Jan  1,  the  mean 
s 

sidereal  time  at  (N^,  T^)  is  given  by 
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T = T (1.0)  + X - 24*^  + 1.0027379093T  + 0 . 0657098222N^ . (8a) 

s s u D 

As  an  example,  for  1977  T (1.0)  = 6^42™7?127  so  if  X = T = 0,  N = 365 
^ s u D 

the  formula  yields  T = 6^41™9?833  while  the  mean  sidereal  time  at  O^UT 

s 

Dec  32,  1977  is  6^41™9f835. 

To  compute  the  fraction  of  the  year  elapsed  to  (from)  the  beginning 
of  the  nearest  Besselian  solar  year  we  proceed  in  two  steps.  Let  x(l.O) 
be  the  fraction  of  the  year  to  Jan  1.0.  Compute  x'  via 

x'  = x(l.O)  + (Nj^  + T^/24)/365.2422.  (9) 

If  N,  < N^(Jul  1),  X = x'.  If  > N,^(Jul  1),  X = X’  - 1.  Thus,  at 
d — D D D 

(334^^,  0*^)  in  1977  x'  = [0.537  + (334  + 0/24)  1/365.2422  = 0.91593  so 
X = -0.08407  since  1977.0  = Jan  0^463. 

Equation  (8a)  is  useful  if  T^(l.O)  is  known  for  any  year.  The  mean 
sidereal  time  at  O^UT  of  any  calendar  date  is  defined  to  be 

T = [6^38"'45!836  + 8640184!542T  + o!o929T^  ]mod  (86400®)  (8b) 

s 

where  T denotes  the  number  of  Julian  centuries  of  36525  days  which,  at 

the  midnight  beginning  the  day,  have  elapsed  since  12^UT  Jan  0,  1900  at 

Greenwich.  Hence,  to  obtain  T^  at  Jan  1.0,  1977  = O^UT  Jan  1,  1977  we 

observe  that  (77  years  of  length  365  days  have  elapsed)  + (19^^  for  leap 

d h h 

years  [1900  was  not  a leap  year])  + (0.5  from  12  UT  to  0 UT  on  Jan  0, 
1900)  = 28124^5.  Straightforward  substitution  into  Eq.  (8b)  yields  T^  = 
6^42™7®127  as  noted  above. 
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Now  that  we  know  the  instant  of  time  that  the  observations  are  to 


be  made  at,  we  can  reduce  the  apparent  place  to  the  topocentric  place. 
For  stars  this  includes  the  effects  of  geocentric  parallax,  diurnal 
aberration,  and  astronomical  refraction.  The  effects  of  geocentric 
parallax  are  less  than  those  of  annual  parallax  by  - (radius  of  earth) /l 
A.U.  Hence,  they  are  universally  ignored.  The  computation  of  diurnal 
aberration  requires  the  observer's  geocentric  coordinates. 

One's  geocentric  coordinates  are  geocentric  distance,  p,  in  units 
of  the  earth's  equatorial  radius  (6378.160km),  geocentric  latitude,  ()>', 
and  geocentric  longitude,  A.  One's  geodetic  coordinates  are  height 
above  mean  sea  level,  H,  in  km,  geodetic  latitude,  (p,  and  geodetic 
longitude,  A.  If  H = 0 then 


psin({)' 

= Ssincf),  pcos(J)'  = Ccos(f), 

(10a) 

tancf)' 

= (1  - f)^tan((). 

(10b) 

A = A, 

(10c) 

P^  = [ 

2 4 2 2 

cos  (})  + (1  - f)  sin  ())]C  , 

(lOd) 

where  f is  the  flattening  of  the  earth  (1/f  = 298.25)  and  the  auxiliary 
quantities  S and  C are  given  by 


1/C^  = cos^cj)  + (1  - f)^sin^4),  (11a) 

S = (1  - f)^C.  (11b) 
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Equations  (11a,  11b,  10b,  and  lOd)  may  be  expressed  directly  in  terms  of 


f.  Through  all  terms  of  the  fourth  order  the  result  is 

C = 1 + f/2  + 5f^/16  + 7f^/32  + 169f^/1024 

- [f/2  + f^/2  + 27f^/64  + llf^/32]cos24) 

+ [3f^/16  + 9f^/32  + 77f^/256]cos4(}) 

- [5f^/64  + 5f'^/32]cos64)  + (35f ^/1024)cos8c() , 

S = 1 - 3f/2  + 5f^/16  + 3f^/32  + 41f'^/1024 

- [f/2  - f^/2  - 5f^/64]cos2<f) 

+ [3f^/16  - 3f^/32  - 19f'^/256]cos44) 

o / 

- (5f  /64)cos6<()  + ( 35 f 71024) cos8(j), 

4)’  = 4)  - [f  + f^/2  - f'^/4]sin24) 

+ [f^/2  + f^/2  + f"^/8]sin44> 

- [f^/3  + f^/2]sin64>  + (f^/4)sin84>, 

p = 1 - f/3  + 5f^/16  + 5f^/32  + 6009f‘^/1024 
+ [f/2  - 13f^/64  - 509f"^/64]cos24) 

- [5f^/16  + 5f^/32  - 517f'^/256]cos44) 

+ [13f^/64  + 13f'^/64]cos64>  - (141f '^/1024)cos84). 


(12a) 


(12b) 


(12c) 


(12d) 


11 


Tf  a correction  for  H is  necessary  then  it  may  be  approximately 
included  by  rewriting  Eq.  (10a)  as 

psin*'  = (S  + H/a  )sinA, 
e 

pcoscf''  = (C  + H/a^)cosi}i, 

where  a is  the  equatorial  radius  of  the  earth, 
e 

The  observer's  speed  relative  to  the  earth's  axis  of  rotation  (due 
to  tlie  rotation  of  the  earth)  is 


(13a) 

(13b) 


V 


(1  revolution/sidereal  day)pa^cos(j>' . 


(14) 


As  there  are  86164.09054  mean  solar  seconds  in  one  mean  sidereal  day 
(see  above),  v = 0. 4651028pcos<J) 'km/sec . Division  of  v by  the  speed  of 
light  in  vacuo  (2.997925  x lO^km/sec)  yields  v ==  0'.'320002pcos<}) ' . The 
corrected  coordinates  (a',  6')  are  given  by 


a'  = a + (vsinl" ) coshsec6 , 
6'  = 6 + (vsinl")slnhsin5 , 


(15a) 

(15b) 


where  h = T - a and  vsinl"  = 1.551416  x 10 
s 

The  last  correction  is  for  astronomical  refraction  and  is  computed 


bv 
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6 , = 6 ' + R' cosn ' , 

obs 

= 6'  + R' sec6 ' cscz' [sincj)  - sin6'cosz']»  (16a) 

a , = a'  + R'sec6  , sinri', 

obs  obs 

= a'  + R'sec6^^^cscz'cos(j)sinh' , (16b) 

where  h'  = - a',  z'  is  the  zenith  distance  and  ri'  is  the  parallactic 

angle  corresponding  to  (a'.  S').  For  any  right  ascension  a,  declination 
6,  and  geodetic  latitude  (j), 


cosz  = sin4>sin6  + cos(j)cos6cosh,  (17a) 

sinzcosp  = sgn(4>)  [sin<j)cos6  - cos4>sin6cosh] , (17b) 

sinzsinri  = cos4>|sinh|,  (17c) 

h = T - a.  (17d) 

s 


The  quantity  R'  (in  radians  above,  seconds  of  arc  below)  is  the 

astronomical  refraction,  R,  corrected  for  the  local  pressure  in  mbar , P, 

and  the  local  temperature  in  °F,  T . 

r 

R'  = (0.5020098P)R/(460  + T^) . (18) 

r 

The  astronomical  refraction  is  the  difference  between  the  true  and 
observed  zenith  distances  (assuming  refraction  to  be  the  sole  cause  of 
the  difference)  and  may  be  approximately  expressed  by 
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R = z' 


R^tanz  , 

1 obs 


(19) 


- z 


obs 


- '*2“"  ^ob3- 


where  R^  and  R^  are  constants.  The  values  used  here  are  R^  = 58V294, 

R„  = 0V0668  and  correspond  to  standard  values  of  (P , T ) = (1015. 92mb, 

2.  r 

50°F).  Equation  (19)  is  accurate  for  z'  £ 75°  and  Includes  effects  due 
to  the  curvature  of  the  earth. 

The  fact  that  only  local  meteorological  conditions  affect  the 
refraction  can  be  simply  deduced  if  the  atmosphere  is  plane-parallel. 
Similarly,  the  leading  term  in  Eq.  (19)  may  be  easily  derived  for  a 
plane-parallel  atmosphere.  These  results  are  also  valid  to  first  order 
since  (height  of  the  atmosphere/radius  of  the  earth)  is  a quantity  of 
the  first  order.  That  the  results  are  valid  to  the  second  order  too  is 
known  as  the  theorem  of  Oriani  and  Laplace. 

We  now  briefly  discuss  the  accuracy  of  these  procedures.  The  mean 
place  to  apparent  place  reduction,  using  J and  J',  and  restricting 
|t|  to  _<  0.5  is  ±0'.'003  in  each  coordinate  (Porter  and  Sadler^).  The 
neglect  of  geocentric  parallax  is  permissable  except  at  the  poles  (e.g., 
within  0?5  of  the  poles  it's  still  less  than  iO'.'Ol) . The  astronomical 
refraction  correction  here  includes  all  terms  through  the  second  order 
and  can  systematically  affect  the  reductions  (because  of  a poor  choice 
of  R'/R,  Rj^,  and  R2)  • However,  a systematic  bias  in  zenith  distance  is 
easily  discovered  and  corrected  for.  Finally,  one  should  use  the  apparent 
sidereal  time  to  compute  the  refraction  correction,  not  the  mean  sidereal 
time.  Except  for  zenith  distances  so  large  that  Eq.  (19)  is  in  doubt, 
this  produces  no  appreciable  error. 
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III.  DIFFERENTIAL  REDUCTIONS 


It  should  be  clear  now  that  a considerable  amount  of  labor  is 
involved  in  reducing  the  position  of  a star.  Even  with  a modern  hand 
calculator,  it  takes  20-30™  to  perform  the  reduction  to  tO'.'Ol.  To 
reduce  N stars,  all  within  a few  degrees  of  each  other,  special  procedures 
have  been  developed  in  order  to  minimize  what  was  once  hand  labor. 

These  methods  have  the  characteristic  that  all  of  the  stars  are  reduced 
relative  to  the  same  point  on  the  celestial  sphere  and  at  the  same 
instant  of  time.  Hence,  these  procedures  are  referred  to  as  differential. 

In  the  age  of  electronic  computers  one  might  question  the  necessity 
or  value  of  analytical  devices  which  shorten  routine  calculations. 

After  all,  the  procedures  of  § II,  applied  to  each  of  the  N stars  individually, 
will  unquestionably  result  in  greater  accuracy  and  easier  programming. 

The  answer  is  partly  psychological,  connected  with  the  esthetics  of  the 
astronomer  versus  those  of  the  arithmetic  registers  of  a central  processing 
unit,  and  partly  practical.  The  practical  side  deals  with  the  ultimate 
accuracy  needed,  the  ability  to  precisely  ascertain  truncation  and 
roundoff  error,  and  the  development  of  widely  applicable  analytical 
tools . 

In  the  remainder  of  this  section  N stars  occupy  a small  5°)  area 

on  the  sky.  They  are  observed  nearly  simultaneously.  When  they  are 

near  the  celestial  poles  additional  labor  is  required  (cf.  § IIIC)  . The 

center  of  the  field  occupied  by  the  stars  has  mean  coordinates  (<a^>, 

<5  >)  and  proper  motions  (<y>,  related  to  those  of  the  stars  via 

o 
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(20a) 


<a  > = Ea  /N,  <6  > = E6  /N, 

o o o o 

<y>  = Ey/N,  <y’>  = Ey'/N. 


(20b) 


The  sums  in  Eqs.  (20)  extend  over  all  N stars  and  parallax  is  neglected. 
The  mean  sidereal  time  of  the  reduction,  T^ , is  the  midpoint  of  the 
observing  interval. 

A.  Star  Constants  and  Besselian  Day  Numbers 

Consider  a position  with  mean  equatorial  coordinates  (<a  >, 

o 

<6  >)  and  a second,  nearby  one,  with  mean  equatorial  coordinates  (a  , 
o o 

6^).  Let  the  associated  proper  motions  be  (<y>,  <y’>)  and  (y,  y'). 

Then  their  apparent  right  ascensions  (declination  is  treated  analogously) 
from  Eq.  (la)  are 


<a>  = <a  > + sinl"[Aa(<a  >,  <5  >)  + Bb(<a  >,  <6  >)  + 
o o o o o 

+ T<y>] , 


(21a) 


a = a + sinl"[Aa(a  , 5 ) + Bb(a  , 6 ) + ... 
o o o o o 

+ Ty], 


(21b) 


where  the  explicit  dependence  on  right  ascension  and  declination  of  the 

star  constants  in  right  ascension  is  indicated.  If  Aa  = a - <a  >,  Aa 

o o o 

= a - <a>,  etc.,  then 


Aa  = Aa  + slnl"{A[a(a  , 5 ) - a(<a  >,  <6  >)]  + 
o o o o o 

+ T[y  - <y>]}. 


(22) 
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However, 


a(a  , 6 ) - a(<a  >,  <6  >)  + y) 

o o o o dx 


Aa  +l5lSi-!i 

o oy 


(<a  >,  <6  >) 
o o 


A6  , (23) 


(<a  >,  <6  >) 
o o 


so 


(24a) 
(24b) 

where  Ay  = y - <y>,  Ay'  = y'  - <y’>,  and 


a - <a>  + Aa  + pAa  + qA6  + rAysinl", 
o o o 

6 = <6>  + A6  + rAa  + sA5  + xAy'sinl", 
o o o 


p = sinl"{ [Acos<a  > - Bsin<a  >ltan<6  > 
o o o 

- [Csin<a  > - Dcos<a  >]sec<6  >} 
o o o 


2 

q = sinl"{ [Asin<a  > + Bcos<a  >]sec  <6  > 
o o o 

+ [Ccos<a  > + Dsin<a  >]sec<6  >tan<6  > 
o o o o 

2 

+ 2Jsec  <6  >tan<6  >}, 
o o 


(25a) 


(25b) 


r = -sinl"{ [Asin<a  > + Bcos<a  >] 


+ [Ccos<a  > + Dsin<a  >]sin<6  >}, 
o o o 


(25c) 


s = -sinl"{Ctanesin<6  > + [Csin<a  > - Dcos<a  >]cos<6  > 

o o o o 

+ 2J'sec^<6  >}. 

o 


(25d) 


17 


Hence,  as  long  as  Aa  , A6  are  not  too  large  (<  0.1  rad  = 5?7) , the  mean 
place  of  any  nearby  star  can  be  reduced  to  its  apparent  place  once  the 
center  of  the  field  has  been  reduced.  As  p,  q,  r,  and  s only  depend  on 
the  position  of  the  center  of  the  field,  it  is  little  more  labor  to 
reduce  all  N stars  simultaneously. 

B.  Independent  Day  Numbers 

The  analytical  complexities  of  the  differential  reduction  can 
be  minimized  if  independent  day  numbers  are  used.  This  is  because, 
unlike  the  star  constants,  the  Independent  day  numbers  do  not  depend  on 
the  position  of  the  center  of  the  field.  Thus,  proceeding  analogously 
to  the  above  [and  using  Eq.  (2a)  in  place  of  Eq.  (la)]  we  have 

<a>  = <a  > + sinl"[f  + gsin(G  + <a  >)tan<6  > + ... 

O L o o o 

+ T<y>] , (26a) 

a = a + sinl"[f  + gsin(G  + a )tan6  + ... 
o o o 

+ xy] , (26b) 


so 


Aa  = Aa  + sinl"{g[sin(G  + a )tan6 
o o o 

+ . . . + T[y  - <y>] } , 


sin(G  + <a  >)tan<6  >] 
o o 

(27) 


or 
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a 


(28a) 


<a>  + Aa  + PAa  + QA6  + TAysinl", 
o o o 

6 - <6>  + A6  + RAa  + SA6  + TAy'slnl",  (28b) 

o o o 

where, 

P = slnl"[gcos(G  + <a  >)tan<6  > + hcos(H  4-  <a  >)sec<6  >],  (29a) 

o o o o 

2 

Q = sinl"[gsln(G  + <a  >)sec  <6  > + hsln(H  + <a  >)sec<6  >tan<6  > 

o o o o o 

+ 2Jtan<5  >sec^<6  >],  (29b) 

o o 

R = -slnl" [gain (G  + <a  >)  + hsln(H  + <a  >)sin<6  >],  (29c) 

o o o 

S = sinl"[hcos(H  + <a  >)cos<6  > - lsin<6  > + J'sec^<6  >].  (29d) 

o o o o 

C.  The  Polar  Regions 

Right  ascension  and  declination  form  an  orthogonal  curvilinear 
coordinate  system,  on  the  celestial  sphere,  with  a non-trivial  metric 
tensor.  Hence,  we  cannot  expect  approximate  procedures  to  work  throughout 
the  right  ascension,  declination  ranges.  However,  (acos6,  6)  do  form  a 
set  of  coordinates  appropriate  everywhere  on  the  celestial  sphere. 

Wherever  cos6<<!  (and  6 = ±81°  is  traditionally  the  dividing  point  in 
astrometry)  it  is  simpler,  and  more  accurate,  to  use  the  direction 
cosines  (cosacos5,  sinacos6,  sin6)  in  place  of  (a,  6).  Thus,  pole  star 
tables  are  frequently  given  in  this  form.  The  above  reduction  can  be 
extended  to  the  direction  cosines  (with  considerable  analytical  complexity) . 
Since  we  never  need  to  observe  artificial  satellites  this  close  to  the 
poles,  the  results  are  not  presented  here. 
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As  a general  rule,  any  approximation  procedure  used  in  spherical 
astronomy  deteriorates  systematically  with  declination.  Thus,  the  position 
dependence  of  the  second  order  terms  in  Eqs.  (1)  is  purely  declination, 
etc.  Hence,  one  must  be  aware  of  systematically  biasing  one's  data 
reduction  in  this  manner. 

D.  Apparent  Place  To  Topocentric  Place 

Once  (<a>,  <6>)  have  been  obtained  they  must  be  corrected  for 
diurnal  aberration: 

<a'>  = <a>  + (vsinl")cos<h>sec<6> , 

<6'>  = <6>  + (vsinl")sin<h>sin<6>, 

<h>  = T - <a>. 
s 

The  (approximate)  differential  reduction  is  obtained  from 

a'  = <a'>  + Aa  + (vsinl")sec<6>[Aasin<h>  + A6cos<h>tan<6>] , (31a) 

5'  =:  <6'>  + A6  + (vsinl")  [-Aacos<h>sin<6>  + A6sin<h>cos<6>] , (31b) 

where,  as  above,  Aa  = a - <a>,  etc.  The  last  correction  is  for  refraction. 


. > 
obs 

= <6'>  + R'cos<r|'>, 

(32a) 

<a  , > 
obs 

= <a'>  + R' sec<6  , >sin<n'>, 
obs 

(32b) 

and,  with  Aa'  = a'  - <a'>,  A6'  =6'  - <6'>, 


(30a) 

(30b) 

(30c) 
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(33a) 


a . ^ <a  . > + Aa'  + lAa ' + JA6 ' , 

obs  obs 

6 ^ =:  <6  ^ > + A6'  + KAa'  + LA6’.  (33b) 

obs  obs 

The  constants  I-L  are  computed  assuming  R = R^tanz^^^  [rather  than  using 
Eq . (19)],  = R'/R,  and  R^  = 58V2  (the  error  is  third-order).  The 

results  are 

I = -R' [1  - tan<z ' >cos<n ' >tan<6 ' > + tan^<z ' >sin^<n ' >] , (34a) 

o 

J = -R^tan<z ' >sin<ri ' >sec<6 ' >[  tan<z' >cos<ri ' > - tan<6’>],  (34b) 

K = -R^tan<z ' >sin<n ' >cos<6 ' >[  tan<z' >cos<f| ’>  + tan<6’>],  (34c) 

and, 

L = -R'[l  + tan^<z'>cos^<n'>] • (34d) 

o 

The  parallactic  angle  and  zenith  distance  are  computed  from  Eqs.  (17). 

Hence,  we  now  have  the  ability  to  completely  and  simultaneously 
reduce  the  positions  of  N stars  (relative  to  any  other  position)  with 
a minimum  of  computation.  We  also  see  that  the  correction  for  astronomical 
refraction  involves  the  most  labor.  It  is  for  this  reason  that,  traditionally, 
this  stage  of  the  reduction  was  left  out.  It's  effects  were  absorbed  in 
the  analysis  of  the  photographic  plate  (see  § V) . In  addition,  because 
ic  is  small,  the  annual  and  diurnal  aberrations  were  not  accounted  for.  The 
complete  reduction  is  recommended. 
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IV.  IDEAL  ASTRONOMICAL  PHOTOGRAPHY 


A.  Standard  Coordinates 

Consider  the  ideal  refracting  telescope  depicted  in  Fig.  1. 

OC  is  the  optical  axis  of  the  telescope  and  GH  is  the  focal  plane.  The 
focal  plane,  normal  to  OC,  contains  a photographic  plate.  OC  intersects 
the  plate  at  its  center.  Produce  OC  to  A where  it  intersects  the  celestial 
sphere.  If  there  were  a star  at  A its  light  would  be  focused  at  0 while 
a nearby  star  at  B (on  the  celestial  sphere)  would  have  its  light  focused 
at  R (on  the  plate).  To  determine  the  relationship  between  the  linear 
size  of  the  photographic  plate,  Z,  and  the  angular  distance  between  the 
corresponding  points  on  the  sky,  L,  Imagine  R to  be  at  the  edge  of  the 
plate  and  the  plane  RCO  to  be  parallel  to  the  edge  of  the  plate.  Then 
plane  trigonometry  in  ARCO  yields 

Z = 2ftanL,  (35) 

where  f is  the  focal  length  of  the  object  glass.  Hence,  from  the  focal 
length  and  the  size  of  the  photographic  plate  one  can  compute  (in  angular 
measure)  the  area  of  the  sky  which  can  be  recorded.  It  is  also  important 
to  know  the  plate  scale  (=  2L/Z)  so  that  a linear  separation  on  the 
plate  can  be  directly  transformed  into  an  angular  separation  on  the  sky. 

Again  referring  to  Fig.  1,  the  plane  tangent  to  the  celestial 
sphere  (also  called  the  plane  of  the  sky)  at  A is  Indicated.  This  plane 
is  (necessarily)  perpendicular  to  OCA  and,  hence,  parallel  to  both  the 
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Fig.  1.  Refracting  telescope  diagram  for  astronomical  photography.  Not 
to  scale. 
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focal  plane  and  the  photographic  plate.  Two  stars  at  I and  B are  shown 
projected  onto  the  tangent  plane  at  K and  D.  Let  <5  OCJ  = 4>  = <5  ACl. 
Then, 


tan({)  = OJ/OC  = AK/CK,  (36) 

so  that  there  is  a similarity  between  the  configurations  of  the  stellar 
images  on  the  plane  of  the  sky  and  on  the  photographic  plate.  If  AS', 

AR'  and  OS,  OR  define  the  positive  directions  for  Cartesian  axes  in  the 
two  planes  (note  the  direction  reversal)  then  coordinates  (?',  o')  on 
the  tangent  plane  are  related  to  coordinates  (^,  o)  on  the  plate  by  the 
scale  factor  AC/OC,  viz. 

= (AC/OC)?,  n'  = (AC/OC)n.  (37) 

The  standard  coordinates,  (?,  n) . (introduced  by  Turner^)  are  measurable. 

.Uist  as  important,  they  can  be  computed  from  the  right  ascension  and 
Jeclination  of  the  point  A (i.e.,  the  point  of  tangency)  and  the  point 
on  the  celestial  sphere,  to  which  they  refer.  I now  demonstrate  these 
statements. 

Figure  2 is  a portion  of  Fig.  1 drawn  from  a different  perspective. 
The  arc  AB  is  that  portion  of  the  great  circle  through  the  point  of 
tangency  and  the  star  at  B.  All  radii  (of  the  celestial  sphere)  connecting 
C to  points  on  AB  lie  in  the  place  of  the  great  circle  through  A and  B 
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Fig.  2.  Same  as  Fig.  1 seen  from  a different  perspective. 
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and  project  onto  the  straight  line  segment  AD  (which  lies  in  the  tangent 
plane) . If  P is  the  North  Celestial  Pole  then  AP  is  the  meridian  for  A. 
Also  drawn  are  n'  axes  with  the  right  ascension  increasing  with 
increasing  Hence,  in  the  figure,  B is  east  of  A's  meridian.  Since 

AD  lies  in  the  tangent  plane  it  is  perpendicular  to  AC.  Similarly,  AU 
is  orthogonal  to  the  great  circle  arc  AP  at  A.  Finally,  the  dihedral 
angle  UAD  is  equal  to  <0  PAB  in  the  spherical  triangle  PAB . Therefore, 
the  projection  (onto  the  plane  of  the  sky)  of  the  arcs  of  great  circles 
preserves  angles.  Let^AB^=  <p,  <J  BAP  = 9 = <^  UAD.  Drop  perpendiculars 
from  D to  AU  (i.e.,  FD)  and  from  D to  AV  (i.e.,  DE) . Then,  from  the 
plane  right  triangles  FDA,  DEA, 

= FD  = ADsinB, 
n'  = ED  = ADcosS. 

l(wever,  AD  = ACtan(|)  so  [from  Eqs.  (36,  37,  38)], 

5 = ftan(J)sin9,  (39a) 

n = ftanc()cos9.  (39b) 

Since  the  focal  length  of  the  object  glass,  f,  merely  serves 
to  define  the  linear  scale,  we  may  set  f = 1 without  loss  of  generality. 
To  complete  the  demonstration  of  the  assertions  made  above,  let  the 
equatorial  coordinates  (relative  to  some  equator  and  equinox)  of  the 


(38a) 

(38b) 
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point  of  tangency  be  (A,  A)  and  those  of  the  star  at  B (relative  to  the 
same  equator  and  equinox)  be  (a,  6) . Then  from  the  spherical  triangle 
PAB  (in  which '"aP^=  90°  - A,'^P'=  90°  - 6,  < APB  = a - A)  , 

cos(()  = sinAsinA  + cosAcosAcos (a  - A),  (40a) 

sin(f)sin6  = cosAsin(a  - A),  (40b) 

sin<})cos6  = sinAcosA  - cosAsinAcos (a  - A).  (40c) 

Hence, 

^ = cotAsin(a  - A)/[sinA  + cotAcosAcos (a  - A)], 

= cosqtan(a  - A)sec(q  - A),  (41a) 

r)  = [cosA  - cotAsinAcos  (a  - A)]/[sinA  + cot  AcosAcos  (a  -A)], 

= tan(q  - A) . (41b) 

Cotq  = cotAcos(a  - A)  and  q is  the  declination  of  that  point  on  the 

celestial  sphere  where  the  great  circle  arc  drawn  from  B intersects  AP 

in  a right  angle.  The  inverse  relationships  are 

tan(a  - A)  = ^secA/[l  - ntanA],  (42a) 

cotAsin(a  - A)  = ^secA/[n  + tanA],  (42b) 

cotAcos(a  - A)  = (1  - r]tanA)/[ri  + tanA].  (42c) 

If  [oi  - a]  is  small  one  should  use  Eq.  (42d)  in  place  of  Eqs . (42b,  c)  . 
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sin6  = [sinA  + r|cosA]/[l  + 


(42d) 


It  has  now  been  demonstrated  that  from  the  (measurable)  standard 


coordinates  plus  the  equatorial  coordinates  of  the  point  of  tangency 
(which  was  chosen  by  the  observer)  one  can  compute  the  right  ascension 
and  declination  of  any  other  point  on  the  photographic  plate.  Conversely, 
from  the  equatorial  coordinates  of  the  star  and  the  point  of  tangency 
one  can  predict  the  standard  coordinates  of  the  corresponding  image 
point  on  the  photographic  plate.  Before  demonstrating  the  utility  (cf. 

§ IVB)  of  this  result  various  series  expansions  derived  from  the  rigorous 
formulas  are  listed  for  reference. 


^ - (a  - A)cosA  - (a  - A) (6  - A)sinA  + (a  - A) ^cosA(3cos^A  - l)/6 


+ . . . , 

(a  - A)cos6  + (a  - A)^cos6(3cos^6  - l)/6 
+ (1/2) (a  - A) (6  - A)^cos6  + ..., 


• • • > 


(43a) 


n - (6  - A)  + (1/4) (a  - A)^sln2A  + (1/2) (a  - A) ^(6  - A)cos2A 
+ (6  - A)^/3  + ..., 

(6  - A)  + (1/4)  (a  - A)^sin26  + (6  - A)^/3  + ..., 


(43b) 


and,  inversely. 


(44a) 
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6 _ A n - (?^/2)tanA  - (C^/2)nsec^A  - n^/3  + ..., 

n - (?^/2)tan6  - n^/3  - ...  . (44b) 

Finally,  note  that  the  convergence  of  Eqs.  (43,  44)  deteriorates  rapidly 
as  |6|  ->■  90°  and  that  |q  - A|  is  always  small.  To  exploit  this  note 

tanq  = tan6sec(a  - A),  (45a) 


and , 


tan(a  - A)  = ?secqcos(q  - A). 


(45b) 


Hence, 


tan(q  - A)  = tan^[(a  - A) /2] sin2q{ [1  - tan^[(a  - A)/2]cos2q},  (45c) 


so 


q A + tan^[(a  - A)/2]sin2q  + (l/2)tan^[(a  - A)/2]sin4q  + ...,  (46) 

which  is  a rapidly  converging  series,  the  truncation  error  being  ^ O'.'Ol  if 
|a  - a|  < 30“. 

B.  The  Method  of  Dependences 

Let  us  now  see  how  the  above  analytical  development  can  be 
used.  Since  only  ideal  astronomical  photography  is  discussed  here. 
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sources  of  error  will  be  glossed  over  (see  § V).  However,  it  should  not 
surprise  the  reader  to  learn  that  the  coordinates  measured  on  the  plate, 
denoted  by  (x,  y) , are  not  exactly  equal  to  the  standard  coordinates. 
Nonetheless,  the  assumption  that  the  relationship  between  their  differences 
is  linear  is  usually  sufficient.  Hence,  we  write  (here) 


^ - X = a^  + bn  + c,  (47a) 

n-y=A^+Bn+C.  (47b) 


The  constants  a,  b,  c.  A,  B,  and  C are  called  plate  constants.  In 
general  they  should  be  (and  will  be)  determined  by  some  estimation 
procedure  such  as  maximum  likelihood  or  least  squares. 

The  problem  we  face  is  the  following:  An  artificial  satellite 

and  N reference  (or  comparison)  stars  have  been  observed  (i.e.,  photographed). 
We  know  the  equatorial  coordinates  of  the  tangential  point  (A,  A),  and 
the  equatorial  coordinates  for  all  of  the  stars,  ^ Hence,  for 

star  j one  can  compute  its  standard  coordinates,  j ^ have  also 

measured  the  coordinates  for  the  stars,  satellite, 

(X,  Y).  Using  a model  for  the  errors  we  wish  to  determine  the  standard 
coordinates  for  the  satellite,  (5,  H)  , and  thence  its  equatorial  coordinates. 

A general  procedure  to  solve  this  problem  was  developed  by 
Schleslnger^ . It  is  called  the  method  of  dependences.  The  minimal  case 
N = 3 illustrates  the  essential  features  of  the  problem  and  can  also  be 
(approximately)  solved  by  graphical  means.  It  also  leads  to  the  optimal 
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configuration  of  the  comparison  stars  relative  to  the  satellite.  For 


these  reasons  we  now  describe  it  in  some  detail. 

The  situation,  on  the  photographic  plate,  is  Illustrated  in 
Fig.  3.  The  reference  stars  are  at  S^,  and  while  the  program 
object  (an  artifical  satellite  here)  is  at  C.  If  one  used  Eqs.  (47)  for 
the  plate  model  and  knew  the  plate  constants,  then  it  would  be  a straightforward 
matter  to  compute  the  satellite's  standard  coordinates  from  its  measured 
coordinates.  To  avoid  actually  calculating  the  plate  constants  let  us 
remember  that  the  least  squares  equations  of  condition  (for  the 
treatment  of  ri  follows  analogously)  are 


- Xj  = aCj  + brij  + c,  J = 1,  2,  3 


(48) 


and,  for  the  artificial  satellite. 


E-X  = aH  + bH  + c. 


(49) 


Multipliers  D.,  j = 1,  2,  3,  called  dependences,  are  introduced  such 


J 


that 


j = 1 


3 

E D (K  - X ) - (H  - X)  = 0. 
_ 1 J J J 


(50a) 


This  yields  three  equations  (equivalent  to  the  normal  equations  of 


least  squares  when  the  coordinate  variances  of  the  reference  stars  are 


equal)  for  the  dependences,  viz. 
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18-9-5719 


Fig.  3.  The  exposed  photographic  plate  with  reference  stars  at  Sj^,  S2 , 
and  and  program  object  at  C. 
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(51a) 


3 


3 

2 D.n.  = H, 

j = 1 1 J 

3 

T.  D.  = 1. 

i = 1 ' 


(51b) 

(51c) 


It  also  follows,  from  the  definition  of  dependences  , that 


3 

T.  D.x.  = X,  (52) 

j = 1 ^ ^ 


and 


3 

c = X + E D. (C.  - X.) . (50b) 

j _ 2^  J J J 

Now,  if  the  linear  plate  model,  Eqs.  (47),  is  sufficient,  then 

since  the  unmodeled  terms  are  of  the  second  order,  the  dependences  will 

be  determined  with  sufficient  accuracy  if  the  measured  coordinates  are 

•k 

substituted  for  the  standard  coordinates  in  Eqs.  (51).  If  this  is  done 
then  the  dependences  are  given  by 


D, 


X X2  X3 
Y X3 

111 


X X3  X3 

Y y,  y3 
111 


X X3  x^ 
Y y,  y^ 
111 


^1  ^2  ^3 


^1  ^2  ^3 


(53) 


*This  departure  from  the  usual  procedures  of  least  squares  analysis  implies  that 
the  original  estimation  problem  was  not  well-posed.  See  § VI  too. 
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The  last  determinant  on  the  right  hand  side  of  Eq.  (53)  is  twice  the 
area  of  the  triangle  formed  by  the  reference  stars.  The  other  determinants 
have  analogous  interpretations  leading  to 


area 

ACS^S^ 

1 

n 

area 

AS^S^S^ 

■ ^"l 

area 

ACS^S^ 

CP2 

area 

AS^S^S^ 

- S2P2 

area 

ACS^S^ 

CP3 

area 

AS^S^S^ 

- S3P3 

(54a) 


(54b) 


(54c) 


Hence,  the  computation  of  the  dependences  has  been  reduced  to  counting 
squares  on  a piece  of  graph  paper.  The  satellite's  standard  coordinates 
are  then  calculated  from  Eqs.  (51a,  b)  and  the  problem  is  solved. 

In  addition  to  being  a neat  labor-saving  trick,  the  geometrical 
significance  of  the  dependences,  coupled  with  the  theorem  that  the 
center  of  mass  of  a plane  triangle  with  a uniform  surface  density  lies 
at  the  intersection  of  its  meridians,  shows  that  the  satellite  should  be 
at  the  "center  of  mass"  of  the  reference  stars  for  maximum  accuracy.  In 
the  general  case  of  N ^ 3 one  can  easily  prove  this  result  too  (Plummer^) 
as  long  as  the  stellar  coordinate  variances  are  equal  (i.e. , the  stars 
have  the  same  "mass") . 

To  summarize,  we  now  know  how  to  go  from  the  photographic 
plate  (charge  coupled  device,  charge  injected  device,  etc.)  to  the  sky 
and  back,  and  how  to  arrange  matters  to  achieve  the  best  possible  internal 
accuracy. 
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V.  REAL  ASTRONOMICAL  PHOTOGRAPHY 


The  idealized  situation  considered  in  § IV  will  not  obtain  in 
practice.  Here  I discuss  the  modeling  of  the  most  common  and  important 
sources  of  error.  These  can  be  divided  into  two  groups;  those  associated 
with  the  measurement  of  the  photographic  place  (centering  error,  rotation, 
non-perpendicularity  of  the  axes)  and  those  associated  with  the  telescope 
(centering  error,  tilt.  Incorrect  focal  length,  radial  distortion, 
decentering  distortion,  and  coma).  In  addition,  the  effects  caused  by 
not  including  astronomical  refraction  and  annual  aberration  in  the 
original  reduction  process  are  considered  too. 

Centering  Error  (Translation) 

This  is  the  error  caused  by  a translation  of  the  photographic  plate 
relative  to  the  measuring  device.  The  differences  between  the  standard 
coordinates  and  the  measured  coordinates  will  be 

C - X = a constant  = h,  (55a) 

ri  - y = a constant  = k.  (55b) 


If  the  photographic  plate  is  rotated  relative  to  the  measuring 
device  by  an  angle  (ijj  > 0 for  a counter  clockwise  rotation)  then. 


Rotation 


C - X = (1  - cosiJj)C  - HsiniJ;, 


(56a) 


r|  - y = ^sinijj  + (1  - cosiJj)r|. 


(56b) 
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Non-Perpendicularity  of  the  Axes 


If  the  X and  y axes  are  not  orthogonal  then. 


C - X = ritanijj, 


(57a) 


n - y = (1  - sec4j)n, 


(57b) 


where  \p  is  the  acute  angle  between  the  p and  y axes  (e.g.,  x = ^ - r)tanl|j,  y 
= r|sec4^) . 

Centering  Error 

If  the  optical  axis  of  the  telescope  (produced)  pierces  the  sky  at 
(A  + 6a,  a + 6A)  instead  of  at  (A,  A)  then  the  standard  coordinates  will 
be  incorrect.  From  Eqs.  (40,  41)  the  leading  terms  are 

^ - X = cosA6A  - nslnAfiA,  (58a) 


Tf  it  can  be  assumed  that  the  net  result  of  all  of  the  other  sources  of 
error  can  be  adequately  described  by  a translation,  a rotation,  and  a 
dilation  (so  that  the  plate  model  is  of  the  form 


ri  - y = 6A  + CslnAAA. 


(58b) 


^ = Ax  + By  + C, 


(59a) 


p = -Bx  + Ay  + D, 


(59b) 


with  |a|>>|b|),  then  one  can  simply  and  rigorously  correct  for  decentering 


via 
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(60a) 


(60b) 


8 


See  Eichorn  . 

Tilt 

If  the  photographic  plate  is  tilted  by  an  angle  o)  relative  to  the 
focal  plane  then  the  relationship  between  the  true  distance  (1)  of  an 
object  from  the  optical  axis  and  its  apparent  (measured)  distance  (£') 
is 

£/£'  = sectiJcos (tjj  - oo)  , (61) 

where  :p  is  the  angular  distance  subtended  by  H at  the  center  of  the 
object  glass,  £ = ftan^J.  If  the  angles  are  small, 


so  this  error  is  of  the  second  order.  In  terms  of  the  differences 
between  the  standard  coordinates  and  the  measured  coordinates. 


£ - £'  - wip^. 


(62) 


(63a) 


(63b) 


where  p and  q are  constants. 
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Focal  Length 


If  the  true  focal  length  is  f (instead  of  f)  then  any  linear 
measure,  on  the  plate  will  be  in  error  by 


This  error  is,  therefore,  linear  in  the  standard  coordinates. 

Radial  Distortion 

A net  radial  distortion  in  the  imaging  process  produces  third  order 
differences  of  the  form 


where  R is  a constant. 

Decentering  Distortion 

When  all  of  the  components  of  the  object  glass  are  not  aligned 
properly  the  resulting  imperfection  is  called  decentering  distortion. 
It  results  in  errors  of  the  second  order  adequately  modeled  by 


m - £/f  = (£/f)(i  - f/f). 


(64) 


^ - X = R^(5^  + n^), 

n - y = Rn(^^  + n^) , 


(65a) 


(65b) 


5 - X = 2[P5^  + Q?n]  + P(5^  + n^), 

n - y = 2[P5n  + Qn^]  + Q(5^  + n^), 


(66a) 


(66b) 


where  P and  Q are  constants. 
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Coma 


Coma  is  a result  of  the  dependence  of  the  focal  length  on  the 
apparent  magnitude  of  the  star.  If  m is  the  apparent  magnitude  then, 


C - X = Sm^,  (67a) 

ri  - y = Smr).  (67b) 


The  constant  S is  strongly  temperature  dependent. 

Differential  Refraction 

When  the  original  reduction  is  only  to  the  apparent  place  and  not 
the  topocentric  place,  the  plate  model  absorbs  the  effects  of  both 
astronomical  refraction  and  diurnal  aberration.  However,  unless  the 
field  is  both  small  in  extent  and  tanz  £ 1,  its  approximate  treatment 
requires  a full  second  order  expression  in  both  coordinates.  When  a 
first  order  treatment  is  adequate  it  can  be  modeled  by 

C - X = R'[(l  + X^)?  + X Y n],  (68a) 

O Z 2 Z 

n - y = R' [X  Y C + (1  + Y^)n],  (68b) 

o z z z 

where  (X  , Y ) are  the  standard  coordinates  of  the  zenith.  R (R')  was 
z z o o 

introduced  in  Eqs.  (34). 

Differential  Aberration 

Since  diurnal  aberration  is  so  small  it  is  (formally)  ignored. 
However,  annual  aberration  may  be  introduced  into  the  plate  model  via 
error  terms  of  the  form 
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C - X = A^, 


(69a) 


n - y = An.  (69b) 

When  this  is  done  and  the  astronomical  refraction  correction  left  out 
too,  only  general  precession  need  be  allowed  for. 

To  summarize,  if  the  full  reduction  to  apparent  place  is  made 
initially,  and  the  field  is  not  too  large,  a linear  model  should  suffice. 
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VI.  THE  TELESCOPE  MODEL 


Now  that  we  understand  the  types  of  errors  to  be  expected  we  turn 

to  the  problem  of  constructing,  in  real  time,  "plate  constants"  for  the 

telescope.  For  the  set  of  reference  stars  we  know  the  computed  equatorial 

topocentric  coordinates,  { (oi  . , 6 . )},  and  the  measured  equatorial 

obs  obs 

topocentric  coordinates,  {(a.^,  5.^,)}.  The  latter  are  subject  to  error 

but  the  former  are  not.  For  the  artificial  satellite  we  know  its  measured 

coordinates,  (A..,,  A^),  and  we  want  to  deduce  (A  , , A , ).  The  simplest 

T T obs  obs 


a 


linear  model,  for  declinations  within  ±81  of  the  equator,  relates  e 
= a , - ot_  (e.g.,  5-x)toe  =a.  - <a  , > (e.g.,  5) » = <5  , 

- <6  , > (e.g.,  Ti).  e = c - <c>,  and  e = m - <m>.  Similarly  for  ej, 

= - 6^.  The  quantities  c and  m refer  to  the  color  and  apparent 

magnitude  of  the  reference  stars  (observed  at  a constant  camera  gain 
setting).  The  average  values  are  denoted  by  <c>,  <m> . If  necessary, 
one  can  obviously  extend  the  model  to  include  more  than  one  color. 
Furthermore,  if  the  camera  gain  is  adjusted  during  each  individual 
observation  such  that  the  size  of  the  disc  visible  to  the  observer  is 
always  the  same  (and  small),  then  the  magnitude  term  (but  not  the  color 
term)  may  be  eliminated. 

Only  actual  observations  and  data  analysis  can  determine  the  relative 
importance  of  the  various  terms  as  a function  of  the  size  of  the  field, 
spread  in  apparent  magnitude,  spread  in  color,  etc.  In  order  to  illustrate 
the  type  of  analysis  necessary  the  linear  model  is  assumed.  Thus, 
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(70a) 


+ 

ve 

a 

+ 

we  ^ 
0 

+ 

xe 

c 

+ 

y"m’ 

"6  = U 

+ 

ve 

a 

+ 

We. 

0 

+ 

Xe 

c 

+ 

Ye  . 
m 

We  obtain  estimates  for  u,  v,  . . . , Y by  the  procedures  of  least  squares 

using  equal  weights.  More  accurate  weighting  cannot  be  used  for  two 

reasons.  First,  the  Smithsonian  Astrophysical  Observatory  Catalog 

contains  incorrectly  computed  coordinate  variances.  Second,  insufficient 

measurements  of  m and  c will  be  obtained  to  accurately  estimate  their 

variances.  After  the  plate  constants  are  obtained  one  computes  (A  , 

obs 

A , ) from 
obs 


+ x(C 


(71a) 


A , - A_  = U + V(A  . - <a  . >)  + W(A  . - <6  , >) 

obs  T obs  obs  obs  obs 


+ X(C  - <c>)  + Y(M  - <m>). 


(71b) 


where  M,  C are  the  apparent  magnitude  and  color  of  the  artificial  satellite. 

Equations  (71),  regarded  as  two,  linear,  simultaneous,  inhomogeneous 

equations  in  the  two  unknowns  A , A , have  a determinant  = 1 + (small 

obs,  obs 

quantity).  Hence,  the  number  of  truly  significant  figures  remaining  at 

this  stage  of  the  calculation  is  critical.  Thus,  one  might  want  to 

introduce  a.^,,  5.^,  on  the  right  hand  sides  of  Eqs.  (70)  initially  via 

e E = a_  - <a  , >,  ->■  E^  = 6_  - <6  , >.  This  would  be  analogous 

a a T obs  0 0 T obs  ^ 
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to  replacing  (^,  n)  by  (x,  y)  in  the  N = 3 version  of  the  method  of 

dependences  [cf.  below  Eq . (50b)].  Then,  the  corresponding  version  of 

Eqs.  (71)  would  still  represent  two,  linear,  simultaneous,  inhomogeneous 

equations  in  the  two  unknowns  A , , A , but  the  determinant  of  the 

obs  obs 

system  would  be  exactly  unity.  Hence,  the  need  for  high  precision 

throughout  the  entire  process  would  be  alleviated.  Unfortunately,  the 

estimation  problem  posed  by  this  substitution  (e.g.,  E^, 

^6  ^6^  requires  _a  priori  knowledge  of  the  variances  and  covariances  of 

e , Or.,  E , and  E^.  The  covariances  can  be  expressed  in  terms  of  the 
a 6 a 0 

variances  (given  a weak  assumption  about  the  perpendicularity  of  the 
telescope's  axes)  but  no  estimate  for  the  variances  is  possible.  Indeed, 
the  telescope  model  is  needed  precisely  because  the  variances  are  unknown. 
Thus,  the  substitution  will  not  be  used. 

There  are  still  three  possible  sources  of  systematic  error.  One  is 
that  the  linear  model  may  not  suffice.  In  order  to  decide  this  point  an 
external  estimate  of  the  ultimate  accuracy  is  needed.  This  is  discussed 
below.  Another  source  of  error  is  due  to  the  fact  that  (typically)  the 
reference  stars  will  be  several  magnitudes  brighter  than  the  artificial 
satellite.  Hence,  if  the  camera  gain  setting  is  kept  constant,  one  will 
be  extrapolating  the  magnitude  term  instead  of  interpolating.  This  can, 
and  should  be,  avoided.  The  last  point  concerns  an  error  introduced  by 
the  reduction  procedure  itself  (depending  on  the  color) . 

The  refractivity  of  air  varies  by  about  two  percent  over  the  wavelength 

o 

range  4000-7000A.  When  the  seeing  is  good  the  stars  are  not  seen  as 
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points  but  as  vertical  spectra  with  the  blue  end  up.  The  effect  is 
small  (see  Fig.  4),  the  angular  distance  between  the  red  and  blue-green 


ends  being  0'.'35,  0'.'60,  1'.'04,  2V24  for  zenith  distances  of  30(15)75°. 
Therefore,  if  all  of  the  reference  stars  are  the  same  color,  the  inclusion 
of  a color  term  in  Eqs.  (70)  will  not  only  incorrectly  model  this  effect, 
but  will  also  systematically  bias  the  results.  Ideally  one  should 
tabulate  R^,  R^,  and  R^  [cf.  Eqs.  (19,  34)]  as  a function  of  c and 
differentially  correct  (from  the  beginning)  for  color.  This  would  be  an 
extremely  difficult  thing  to  do.  Nonetheless,  one  can  estimate  this 
error  (see  below) . 

One  can  obtain  an  estimate  for  the  internal  error  in  the  artificial 

satellite's  topocentric  position  as  follows:  Solve  Eqs.  (71)  analytically 

for  A , , A , , say  A , = f(u,  v, . . .,  Y).  The  variance  of  A , is 
obs  obs  obs  obs 

given  by 


where  E is  the  covariance  of  u,  v,  etc.  The  solution  of  the  normal 
uv 

equations  provides,  coupled  with  the  actual  unbiased  residual,  estimates 

for  the  variances  and  covariances  appearing  in  Eq.  (72).  Thus,  an 

estimate  for  the  internal  variance  of  A , (the  situation  is  similar  for 

obs 

the  declination)  is  available.  As  a formal  result  this  is  fine,  but 


obs 


3f  Bf  , . „ 3f  9f 

^ uv  3u  3v  xy  3x  3Y 


(72) 
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Fig.  4.  Dispersion  caused  by  atmosphere  in  the  UBVRI  system. 
Taylor . ^ 


After 
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probably  useless.  It  is  the  external  error  (e.g.,  relative  to  the 
reference  system  of  the  FK4)  that  is  of  interest.  A better  way  to 
estimate  this  is  to  have  observed  and  reduced  N + 1 stars,  allowing  N to 
contribute  to  the  plate  model.  The  extra  star,  which  most  closely 
matches  the  artificial  satellite  in  color,  apparent  magnitude,  and 
position,  serves  as  a control  on  the  entire  procedure. 
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VII.  ARTIFICIAL  SATELLITE  REDUCTIONS 


At  this  stage  of  the  analysis  I assume  that  the  topocentric  equatorial 

coordinates  of  the  artificial  satellite,  (A  , , A , ) are  known  and  one 

obs  obs 

wants  to  compute  geocentric  equatorial  coordinates  (A,  A).  To  do  this 
one  needs  to  correct  for  refraction,  planetary  aberration,  diurnal 
aberration,  and  geocentric  parallax.  Three  of  these  involve  the  distance 
to  the  artificial  satellite.  If  this  is  not  known  it  must  be  estimated. 

An  estimation  procedure  is  described  in  § VIIA  and  the  rest  of  § VII 
discusses  the  four  corrections. 

A.  Distance  Estimation 

In  this  subsection  only  lower  case  variables  refer  to  geocentric 
values  (a,  6,  h = hour  angle,  w = angular  speed,  d = distance)  and  upper 
case  variables  refer  to  topocentric  values.  Also,  in  this  subsection 
only,  the  difference  between  the  two  sets  is  due  solely  to  geocentric 
parallax. 

The  geocentric  location  of  the  observer  is  given  by  ^ = (s  = 

pa  , d)',  T ).  The  fundamental  relationship  between  the  geocentric 
e s 

location  and  the  topocentric  location  is 

d = D + s . (73) 

In  component  form. 
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dcos6cosa  = DcosAcosA  + scoscb'cosT  , (74a) 

s 

dcos6sina  = DcosAsinA  + scoscb'slnT  , (74b) 

s 

dsin6  = DsinA  + ssincf)'.  (74c) 

These  equations  may  be  manipulated  to  yield 

sin(a  - A)  = (s/d) coscf) ' secAsinH,  (75a) 

sin(6  - A)  = (s/d)sln({)  cscTsln(r  - A),  (75b) 

tanP  = tan4/cos[(A  - a)/2]sec[H  + (A  - a)/2],  (75c) 

d = Dsin(A  - r)csc(6  - r),  (75d) 

2 2 2 

d = D + s + 2Ds  (sin({)'sinA  + cos(})'cosAcosH)  , (75c) 

or,  inversely, 

tan(a  - A)  = asinh/(l  - acosh) , (76a) 

tan(6  - A)  = bsin(Y  - <S)/[1  - bcos(y  - 6)],  (76b) 

a = (s/d)coscl)'sec6,  (76c) 

b = (s/d)sin(l) ' cscY,  (76d) 

tanY  = tanF  = tan(J)' cos  [ (a  - A)/2]sec[h  + (a  - A)/2],  (76e) 

D = dsin(6  - y)csc(A  - Y) > (76f) 

2 2 2 

D = d + s - 2ds(sin({)'sin6  + cos(J) ' cos6cosh)  . (76g) 

These  relationships  should  be  used  rather  than  (74)  since  they  will 
always  yield  more  accurate  results. 
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In  order  to  estimate  a distance  one  needs  three  positions  or 
two  positions  and  one  angular  velocity.  Utilizing  the  first  data  set  is 
a standard  problem  in  celestial  mechanics  and  I shall  not  treat  it  here. 
Instead,  I will  develop  the  second  approach. 

Expand  a - A,  5 - A in  powers  of  s/d.  The  first-order  result 
is 


a - A = (s/d)cos(j)'sinHsecA,  (77a) 

6 - A = (s/d)  [sin(j)'cosA  - cos(j>'sinAcosH] . (77b) 

Differentiate  Eqs.  (77)  with  respect  to  time  (rigorously  Ephemeris  time 

2 *2  2 *2  • 

but  mean  solar  time  is  sufficient)  and  form  w = a cos  6+6  where  z = 
dz/dt  V z.  The  leading  terms  are 

00^  = + 2(s/d)  {-fJ^[sin(f)*sinA  + cos(}>'cosAcosH] 

• • • • • 

+ [AcosAcosH  + AsinAsinH]T  cost})'  - (d/d)  [A(sinc|>'cosA 

s 

- cos(() ' sinAcosH)  + Acos4> ' cosAsinH]  } , (78) 

• 2 

so  that  not  only  d but  d/d  are  required  to  estimate  o)  . It  turns  out 
that  d/d  can  be  approximately  computed  without  knowledge  of  d.  I now 
show  this. 

Differentiate  Eqs.  (75a,  75b)  with  respect  to  time.  The  exact 

result  is 
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(79a) 


• • • • • 

d/d  + acot(a  - A)  - 5tan5  = Acot(a  - A)  + HcotH, 

d/d  + otEF  + 6cot(6  - A)  = AEF  + A[cot(6  - A)  - cot(r  - A)] 

+ fetan[H  - (a  - A)/2],  (79b) 


where 


E = sinPcosP [cot (F  - A)  - cotP] , (80a) 

2F  = tan[ (a  - A)/2]  + tan[H  - (a  - A)/2].  (80b) 


In  addition,  by  differentiating  the  geocentric  relationship  which  expresses 
the  constancy  of  the  direction  of  the  total  orbital  angular  momentum 
vector,  one  finds. 


asln5cos6cot (a  - Q')  -5=0, 


(79c) 


where  0,'  is  the  longitude  of  the  ascending  node.  The  only  approximation 
involved  in  determining  d/d  from  Eqs.  (79)  concerns  the  computation  of 
. One  calculates  from  the  two  topocentric  (instead  of  geocentric) 
positions,  viz. 


tanSl' 


tanA^sinA^  - tanA2SlnAj^ 
tanA^cosA2  - tanA^cosA^^ 


(80c) 


Thus,  the  observational  data,  the  geometry,  and  the  physics  provides 
three,  linear,  simultaneous,  inhomogeneous  equations  in  the  three  unknowns 


50 


• * • * 

a,  6,  d/d.  We  now  continue  the  estimation  problem  for  d regarding  d/d 

as  known . 

If  the  orbit  were  circular  and  if  to  = then  d would  be  given 
by 


d = d^  = (81) 


where  G is  the  gravitational  constant  and  M is  the  mass  of  the  earth; 

14  3 2 

GM  = 3.98603  x 10  m /sec  . A better  guess  for  d corrects  d for  the 

w 1 

effects  of  an  eccentric  orbit.  I write  this  as 


= F(e)d^. 


(82) 


F(e)  is  discussed  below.  A new  estimate  for  the  distance,  d^,  is  then 

2 

calculated  using  d = d2  in  Eq.  (78),  and  then  putting  (O  into  Eq.  (81) 

2 

in  place  of  . Finally,  the  last  estimate  for  d is  computed  from  Eq. 

(82)  with  d^  replacing  d^. 

F(e)  is  the  average,  over  an  orbital  period,  of  the  ratio  of 
the  geocentric  distance  in  the  elliptic  orbit  to  the  geocentric  distance 
in  a (fictitious)  circular  orbit  with  the  same  instantaneous  angular 
speed.  Using  standard  relationships  F(e)  can  be  written  as  (e  = eccentricity, 
E = eccentric  anomaly) 


F(e)  = 1(1  - e^)^^^/2TT] 


X 


2tt 

(1  + ecosE)^^^dE 


= P2/3(e)/e^^^;  = 1/(1  - e^) , 


(83) 
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where  P 


is  the  Legendre  function  of  order  2/3.  Table  1 lists  F(e) 


2/3 


for  e = 0(0.05)1.  If  the  eccentricity  can  not  be  estimated  to  ±0.1 
(from  other  information)  then  an  average  value  of  F(e)  should  be  used. 


determined  from  Eqs.  (79)  because  of  the  dominance  of  the  cotangents  in 
Eqs.  (79a,  79b).  Since  it  is  small,  it  may  be  better  to  set  it  equal  to 
zero  in  Eq.  (78)  if  one  doubts  its  sign.  The  overall  accuracy  is  then 
±10  percent  in  d for  d ^ 3s,  increasing  with  d. 

B,  Refraction 

The  light  rays  from  an  extra-terrestrial  object  must  pass 
through  both  space  and  the  earth's  atmosphere  before  reaching  an  earthbound 
telescope.  I assume  that  the  extra-atmospheric  path  is  a straight  line. 
UTien  the  light  enters  the  atmosphere  the  path  becomes  curved,  with  a 
continually  changing  curvature.  The  light  ray's  path  is  bent  towards 
the  geocentric  radius  vector  to  the  source  so  that  refraction  can  alter 
zenith  distance  but  not  azimuth.  The  difference  in  direction  between 
the  final  tangent  vector  to  the  light  ray's  path  (e.g.,  at  the  telescope) 
and  the  initial  tangent  vector  to  the  light  ray's  path  (e.g.,  at  the  top 
of  the  atmosphere)  is  called  the  astronomical  refraction.  In  Fig.  5 
<5  S'OZ  (the  telescope  is  at  0,  Z denotes  the  zenith  of  0)  is  the  apparent 
zenith  distance  of  the  source  and  ^ SO'Z  is  the  true  zenith  distance  of 


o 


(84) 


Experience  with  this  procedure  indicates  that  d/d  is  poorly 
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TABLE  I. 


ECCENTRICITY  CORRECTION  FACTOR 


e 

F(e) 

0 

1 

0.05 

0.999027 

0.10 

0.996101 

0.15 

0.991197 

0.20 

0.984275 

0.25 

0.975276 

0.30 

0.964119 

0.35 

0.950697 

0.40 

0.934873 

0.45 

0.916469 

0.50 

0.895257 

0.55 

0.870938 

0.60 

0.843115 

0.65 

0.811248 

0.70 

0.774575 

0.75 

0.731972 

0.80 

0.681663 

0.85 

0.620560 

0.90 

0.542382 

0.95 

0.429811 

1 

0 
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18-9-5721 


Fig.  5.  Diagram  for  refraction  calculation.  C is  the  center  of  a spherical 
earth,  the  observer  is  at  0,  and  light  from  S enters  the  top  of  the  atmosphere 
at  A.  Not  to  scale. 
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★ 

the  source.  The  difference  between  these  angles  is  the  parallax  of  the 

source  as  seen  from  O'  relative  to  0.  To  determine  the  height  of  O' 

(00'=  L),  regard  the  earth  as  a perfect  sphere  with  radius  a . Then, 

e 

for  a spherical  atmosphere, 


1 + L/a  = u sinz  cscz  = u sin(z 

e o app  true  o true 


R) cscz 


true  ’ 


(85) 


where  is  the  index  of  refraction  of  the  atmosphere  at  the  telescope 

and  R = z - z is  called  the  astronomical  refraction.  The  net 

true  app 

refraction  is  R - r where  r (the  parallactic  refraction)  measures  the 
angle  subtended  by  L along  the  original  light  path  (e.g.,  O'S"  = OS"  in 
Fig.  5).  The  parallactic  refraction  for  the  natural  celestial  objects 
can  only  be  appreciable  for  the  Moon  (and  then  only  at  large  zenith 
distances).  Hence,  it  is  universally  ignored  except  for  the  precise 
predictions  of  solar  eclipses.  However,  r increases  as  the  geocentric 
distance  decreases,  reaches  one  minute  of  arc  for  an  object  500km  above 
the  earth's  surface.  Hence,  its  universal  neglect  for  artificial  satellites 
will  result  in  a systematic  bias. 

An  approximate  expression  for  r (D  = topocentric  distance  in 
km,  H = height  above  earth's  surface  in  km)  is. 


rcotz  cosz  = -(482'.'9/D)  [ 1 - 0.00130(2sec^z 

app  app  app 

+ tan^z  ][1  - exp (-0 . 1205H) ] 
app 

- 0'.'0757 (sec^z  + tan^z  )Hexp (-0 . 1205H)  . (86) 

app  app 


*In  this  subsection  only  the  difference  between  true  and  apparent  angles 
is  due  solely  to  refraction. 
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Thus,  the  equations  analogous  to  Eqs.  (16)  for  artificial  satellites  are 


6 

app 


a 

app 


= 6 + (R'  - r' )sec6  cscz  [sin(j> 

true  true  true 


- sin6  cosz  ], 
true  true 

= a + (R'  - r')sec6  cscz  cosd)sinh 

true  app  true  true 


(87a) 

(87b) 


These  equations  must  be  solved  iteratively  for  (a^  , 6^  ) and  r'/r  = 

^ true  true 

R'  /R. 

C.  Planetary  Aberration 

Planetary  aberration  refers  to  the  aberration  of  light  due 
explicitly  to  the  motion  of  the  source.  Thus,  since  the  artificial 
satellite  is  D km  away,  it  required  D/c  (c  = speed  of  light  = 2.997925  x 
lO^km/sec)  seconds  for  the  light  to  reach  the  observer.  Neglecting  this 
correction  ( - 2"  for  a geostationary  satellite)  will  also  systematically 
bi  is  the  results.  If  (A'",  A"')  are  the  satellite's  topocentric  coordinates 
cc  "rected  for  refraction  and  (A",  A")  are  its  topocentric  coordinates 
corrected  for  both  refraction  and  planetary  aberration  then,  (approximately) 


A"  = A'”  - AD/c, 
A"  = A'"  - AD/c. 


(88a) 

(88b) 


D.  Diurnal  Aberration 

This  is  exactly  as  for  the  stars.  Hence,  (A',  A')  is  related 
to  (A",  A")  via  [cf.  Eqs.  (15)], 
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A"  = 

A' 

+ (vsinl")cosH' 

' secA' , 

(89a) 

A"  = 

A' 

+ (vsinl")sinH' 

sinA' 

(89b) 

H'  = 

T 

s 

- A'. 

(89c) 

E.  Geocentric  Parallax 

The  last  correction  one  must  apply  to  (A',  A')  to  yield  geocentric 
coordinates,  (A,  A),  is  the  one  for  geocentric  parallax.  Adapting  Eqs. 

(75)  to  the  notation  here  we  calculate  successively 

d^  = + s^  + 2D  s(sin(J) 'slnA'  + cos(j)' cosA' cosH ' ) , (90a) 

A = A'  + arcsin[  (s/d)tos({>'secA'sinH' ] , (90b) 

A = A'  + arcsin[  (s/d)sin<{>'cscrsin(r'  - A)],  (90c) 

with 

tanP'  = tancj) ' cos [ (A'  - A)/2]sec[H'  + (A’  - A)/2].  (90d) 
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